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ABSTRACT 

The red spectral shape of the visible to near infrared reflectance spectrum of the sharply-edged ring-like disk 
around the young main sequence star HR 4796A was recently interpreted as the presence of tholin-like complex 
organic materials which are seen in the atmosphere and surface of Titan and the surfaces of icy bodies in the 
solar system. However, we show in this Letter that porous grains comprised of common cosmic dust species 
(amorphous silicate, amorphous carbon, and water ice) also closely reproduce the observed reflectance spectrum, 
suggesting that the presence of complex organic materials in the HR4796A disk is still not definitive. 

Subject headings: circumstellar matter — dust, extinction — infrared: stars — planetary systems: protoplanetary 
disks — stars: individual (HR 4796 A) 



\. INTRODUCTION 

HR4796A is a nearby (distance to the Earth d (fl ± 3pc) 
young main-sequence (MS) star (age w 8 ± 3Myr; Stauffer et 
al. 1995) of spectral type AOV (effective temperature T^a « 
9500 K) with a large infrared (IR) excess which has recently 
aroused considerable interest. Imaging observations describe 
dust scattering in the optical (Debes et al. 2008) and near-IR 
(Augereau et al. 1999; Schneider et al. 1999) and thermal emis- 
sion at mid-IR (Jayawardhana et al. 1998; Koerner et al. 1998; 
Telesco et al. 2000; Wahhaj et al. 2005). They reveal a ring-like 
disk with maximum at ^70AU distance from the central star 
and ^ 17 AU width that is sharply truncated at the inner and the 
outer edge. The structure of the HR4796A disk has important 
implications for planetesimal evolution (Kenyon et al. 1999). 
Furthermore, possibly existing planets may generate disk asym- 
metries through gravitational confinement or perturbation (Wy- 
att et al. 1999) or form rings with sharp edges (Augereau et al. 
1999, Klahr & Lin 2000, Thebault & Wu 2008). 

Very recently, Debes et al. (2008) measured a visible to near- 
IR photometric reflectance spectrum of the dust ring around 
HR4796A. To fit the observed spectrum (which is character- 
ized by a steep red slope increasing from A « 0.5 /im to 1.6 /im 
followed by a flattening of the spectrum at A > 1.6/im), Debes 
et al. (2008) argued for the presence of tholin-like organic ma- 
terial in the disk around HR4796A. Tholin, a complex organic 
material, was detected as a major constituent of the atmosphere 
and surface of Titan and the surfaces of icy bodies in the solar 
system. The detection of tholin in the HR4796A disk - if con- 
firmed - would imply that these potential basic building blocks 
of life may be common in extra-solar planetary systems as well 
(see van Dishoeck 2008 for an overview of organic matter in 
space). Its young age places HR4796A at a transitional stage 
between massive gaseous protostellar disks around young pre- 
MS T-Tauri and Herbig Ae/Be stars (~ 1 Myr) and evolved and 
tenuous debris disks around MS "Vega-type" stars 100 Myr; 
see Jura et al. 1993, Chen & Kamp 2004). Detecting organic 
matter in this system will provide valuable information about 
the formation and evolution of planetary systems. 



The observed thermal emission of the disk, however, was 
previously reproduced with a model population of porous 
grains consisting of the common cosmic dust species amor- 
phous silicate, amorphous carbon, and water ice (Li & Lunine 
2003a). In this Letter we question the existence of tholin as a 
major dust component in the HR4796A disk and study alter- 
native dust models to reproduce the observed reflectance spec- 
trum. 

2. MODEL 

Dust in the disk around HR4796A scatters starlight at visi- 
ble to near-IR wavelengths and emits thermally in the IR. The 
total flux of the disk Fdisk(A) is the sum of the scattered light 
■Psca(A) and the dust thermal emission Ftherm(A). At A < 2.2 fim 
the dominant contribution to Fdisk(A) comes from the starlight 
scattered by dust Fsca(A) which is calculated from 



^:sca(A) = 



— I a{r)2iTrdr 
2r ' ' 



Csca( A, fl) $( A, a) n(a) da 



(1) 



where li « 67 pc is the distance from the star to the Earth; is 
the stellar atmospheric flux approximated by the Kurucz model 
for AOV stars with T^ff = 9500 K, lgg = 4.5 and a solar metallic- 
ity (Kurucz 1979); = 1.7 Rq is the stellar radius; ri„ and rout 
are respectively the inner and outer boundaries of the disk; a(r) 
is the dust surface density distribution; Csca(A,a) is the scatter- 
ing cross section of spherical dust of radius a at wavelength 
A; $(A,fl) is the phase function approximated by the Henyey- 
Greenstein function at a scattering angle of 6 = 90° (see Debes 
et al. 2008) with the asymmetry parameter g calculated from 
Mie theory for dust of size a is the dust size distribution 

which is taken to be a power-law n(a) oc a~" with a lower-cutoff 
flmin, an upper-cutoff Omax, and a power-law index a. We take 
flmax = 1000 /im and treat fl,nin and a as free parameters. 

Foflowing Kenyon et al. (1999), Klahr & Lin (2000), and Li 
& Lunine (2003a), we approximate the dust surface density dis- 
tribution by a Gaussian-type function a{r) = crpexp[-41n2{(r- 
rp)/A}^] where rp is the radial position where a(r) peaks, A is 
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' The asymmetry parameters g of the best-fit models averaged over the size distributions show a gradual decrease from g fa 0.96 at A = 0.55 fj,m to g 0.85 at 
A = 2.5 /im, indicating that the dust is highly forward-throwing. We should note that one should not compare the g values calculated here with that of Debes et al. 
(2008), derived from the observed surface brightness of the disk together with the Henyey-Greenstein phase function. The former specify the degree of scattering in 
the forward direction (9 = 0° ) of the dust (see Li 2008), while the latter approximate the mean scattering properties of the dust averaged over the entire disk. 
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the full width half maximum (FWHM) of the distribution, and 
<Tp is the mid-plane surface density at r = rp. We take rp = 70 AU 
and A = 15 AU (see §2.1 of Li & Lunine 2003a). We take 
rin = 40 AU and rout = lOOAU since there is little dust even at 
r < 55 AU or r > 85 AU (see Fig. 5 of Li & Lunine 2003a). 

We consider porous grains composed of common dust 
species (amorphous silicate, amorphous carbonfl and water ice; 
see Li & Lunine 2003a). The optical properties of porous 
dust are determined with Mie theory in combination with the 
Bruggeman effective medium theory (Bohren & Huffman 1983; 
see eqs.7-9 of Li & Lunine [2003b] for a detailed descrip- 
tion). We take the mass ratio of amorphous carbon to amor- 
phous silicate to be mcai-b/'Wsii = 0.7 and the mass ratio of 
water ice to amorphous carbon and amorphous silicate to be 
wJice/wcarb + 'Wsii = 0.8, as inferred from the cosmic abundance 
constraints (see Appendix A of Li & Lunine 2003a). Porous 
dust models consisting of fluffy aggregates of amorphous sil- 
icate, amorphous carbon and water ice with such mixing ra- 
tios have been shown successful in reproducing the IR to sub- 
millimeter dust emission spectral energy distribution of the 
HR4796A disk (Li & Lunine 2003a, Sheret et al. 2004). 

For amorphous silicate dust, we assume amorphous 
(Mg,Fe)Si04, an amorphous material with olivine-normative 
composition for which we take optical constants from Jager et 
al. (1994). The optical constants of amorphous carbon and wa- 
ter ice are taken from Rouleau & Martin (1991; "AC" type) 
and Warren (1984), respectively. Crystalline silicates are not 
included in our model calculations since observations at 8- 
13/im with the mid-IR Keck LWS show no crystalline sili- 
cate emission features (Kessler-Silacci et al. 2005). Although 
the optical properties of crystalline silicates differ from that of 
amorphous silicates^ the resulting reflectivity is similar. More- 
over, the amount of crystalline silicates, if observed in disks, is 
significantly smaller than that of the amorphous silicates. 

A major characteristic of porous dust is its porosity P (i.e. 
the fractional volume of vacuum in a porous grain). Li & Lu- 
nine (2003a) have shown that the IR emission of the HR 4796A 
disk is best fit by dust with P = 0.90. We should note that the 
best-fit porosity of P = 0.90 refers to the porous aggregates of 
amorphous silicate and amorphous carbon; the porosity is re- 
duced to P Ri 0.73 when ice fills in some of the vacuum un- 
der the assumption of a complete condensation of all condens- 
able volatile elements as ice (see Appendix B of Li & Lunine 
2003a). It is expected that the dust in the HR4796A disk will 
be coated by ice since at ^70 AU from the star the dust will be 
cooler than 110-120K and ice condensation will occur (see 
Li & Lunine 2003a). Assuming a complete ice condensation, 
the original porosity P for the porous mixture of silicate and 
carbon will be reduced to max{0, [1 -2.68(1 -P)]}. Unless 
stated otherwise, in the following the porosity P refers to the 
fractional volume of vacuum in ice-coated porous dust. We will 
consider a range of porosities: P = 0,0.2,0.5,0.73,0.9 (which 
correspond to P w 0.63 , 0.70, 0.8 1 , 0.90, 0.96 if ice is removed). 

3. RESULTS 

We compare Ff^^aW/F^ with the visible to near-IR re- 
flectance spectrum of the HR4796A disk compiled by Debes 

^ The other carbon dust species widely considered in astrophysical modeling are graphite, hydrogenated amorphous carbon, quenched carbonaceous composite, and 
organic refractory. Their optical properties are not qualitatively different from that of amorphous carbon. 

^ While crystalline silicates in astrophysical regions are often found to be Mg-rich and Fe-poor (e.g. see Molster & Kemper 2005), astronomical amorphous silicates 
appear to have a similar fraction of Mg and Fe as implied by the strong UV/visible absorptivity required to model the circumstellar emission (Jones & Merrill 1976, 
Rogers et al. 1983) and interstellar extinction (Draine & Lee 1984). This suggests that amorphous silicates could be less effective in scattering the visible to near-IR 
starhght. 



et al. (2008) either from archival images (HST/STIS 50CCD 
at A = 0.585 Aim, HST/NICMOS FllOW at A= 1.1/im, F160W 
at A = 1.6 i-im) or from their newly obtained HST/NICMOS im- 
ages (F171M, F180M, F204M and F222M at A= 1.71, 1.80, 
2.04, 2.22 jim, respectively). 

By varying Omin and a, we try to fit the observational data 
of Debes et al. (2008) with the porous dust model consisting 
of common dust species for which the mass mixing ratios are 
taken from the cosmic abundance considerations (see §2). We 
first take the porosity to be f = 0.73 (which corresponds to 
P = 0.90 for the porous dust without ice coating; the model 
with P = 0.90 gives the best fit to the observed IR emission, as 
shown in Li & Lunine 2003a). As illustrated in Figure [T] the 
model with P = 0.73, flmin ~ 2 /im, and a « 2.9 closely repro- 
duces the reflectance spectrum of HR4796A. In Figure |2] we 
plot the deviations of the model predictions from the HST pho- 
tometry at each of the seven bands. The deviations are smaller 
than ^7.5% for all wavebands. It is even more encouraging 
that the same model also closely fits the IR emission of this 
disk (see Fig. [3]), except its slight deficiency at ~10/im (but 
still well within the observational uncertainties). The model re- 
quires a total dust mass of ^4.27 x lO^^g and results in a verti- 
cal visible optical depth of ry « 0.048 at r = rp and a mid-plane 
radial optical depth of Ty ~ 0.56. This justifies the optical thin 
approximation employed in this work. 

In contrast, the model with P = 0.73, Omin ~ 1 /^m, and 
a « 2.9 provides the best fit to the observed IR emission (see 
Li & Lunine 2003a and Fig. [3]). However, this model does not 
closely fit the observed reflectance spectrum (see Figs. 11121) . al- 
though the model reflectance spectrum does exhibit a general 
trend (i.e. a steep red slope at ~ 0.5-1 .6 /im and subsequent flat- 
tening off; see Fig.[T]) similar to that observed in the HR4796A 
disk. 

The best fit to the observed reflectance spectrum is given 
by the model with P = 0.50, Omin ~ 1 /^m, and a « 2.8 (see 
Fig.[T]). The deviations from the reflectance data are smaller 
than ^6.2% for all wavebands (see Fig.|2|. However, the model 
does not fit the observed IR emission (see Fig.[3]l. 

Generally speaking, highly porous dust (with P > 0.85) does 
not fit the observed reflectance spectrum well, while more com- 
pact dust (with P < 0.60) is too cold to reproduce the observed 
IR emission. In Table [T] we list the parameters used in the 
different models: the minimal grain radius a^in, the exponent 
of the size distribution a, the porosity P, and x^/dof, where 
dof = A^data-A'para IS the "degree of freedom" (A^data = 7 is the 
number of data points at A = 0.585, I.I, 1.6, 1.71, 1.80, 2.04, 
2.22 /im, A'para = 3 is the number of free parameters: flmin, a and 

py 

4. DISCUSSION 

Debes et al. (2008) calculated the reflectance spectra of "as- 
tronomical silicates", water ice, hematite FeaOs (which is found 
on the surface of Mars and responsible for its redness), and UV 
laser ablated olivine (which has been used to explain the spec- 
tral reddening of silicate-rich asteroids due to space weathering; 
Brunetto et al. 2007). But the scattering spectra of these miner- 
als and ice are too neutral at A ^ 0.5-1.6 /im to match the steep 
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red spectral slope of the reflectance spectrum of the HR4796A 
disk. Debes et al. (2008) therefore resorted to organic mate- 
rials. They found that tholin or its mixture with other dust 
species (e.g. water ice or olivine) are able to reproduce the 
observed red spectral slope. This led them to suggest that ''the 
presence of organic material is the most plausible explanation 
for the observations" , with a cautionary note that "... longer 
wavelength scattered light observations will further constrain 
the ( tholin-based) grain models, particularly around 3.8-4 fim, 
where a large absorption feature is seen for different grain sizes 
of tholins. This would help to directly confirm whether Titan 
tholins are an adequate proxy for the material in orbit around 
HR4796A:' 

However, as shown in Figure [1] simple porous dust models 
consisting of dust species (amorphous silicate, amorphous car- 
bon, water ice) which are commonly considered to dominate in 
the interstellar medium (ISM), envelopes around evolved stars, 
and dust disks around young stars closely reproduce the ob- 
served reflectance spectrum of the HR4796A disk. Our model 
provides at least a viable alternative to the tholin-based models 
of Debes et al. (2008). While the tholin organic dust model pre- 
dicts a strong feature around 3.^-4- fim characteristic of tholin 
(Debes et al. 2008), the porous dust model presented here pre- 
dicts a strong band at ~ 3. 1 /im, attributed to the O-H stretching 
mode of water ice. 

The tholin organics of which the optical constants were 
adopted by Debes et al. (2008) were made from DC discharge 
of 90% N2 and 10% CH4 gas mixture (Khai-e et al. 1984). They 
are extremely N-rich and optically very different from amor- 
phous carbon. The dust in the HR4796A disk should be con- 
tinuously replenished. This is indicated by the low size cutoff 
ctmin (of a fcw micromctcrs) of the dust required to reproduce 
the reflectance spectrum (see §2 and Table [1) combined with 
considerations of dust lifetimes based on the radiation pres- 
sure and Poynting-Robertson effects (see Fig. 9 of Li & Lunine 
2003a). The replenishing source would likely arise from col- 
lisional cascades of larger bodies like planetesimals, asteroid- 
like and comet-like bodies. We suggest, with interstellar dust 
as the building blocks of the parent bodies, it is more reason- 
able to assume that the dust in the HR 4796A disk is cornposed 
of amorphous silicate, amorphous carbon, and water iceQ 

The model which best fits both the observed reflectance spec- 
trum and the IR emission requires highly porous dust (with 
P = 0.73 which corresponds to P = 0.90 if ice is sublimated; 
see §2). While it is natural to recognize that cold conglomera- 



tion of dust grains in molecular clouds can lead to highly porous 
dust structures!] at a first glance, it is harder to accept that com- 
paratively more violent collisions between larger bodies would 
result in such a morphology in the resulting debris. 

To address this concern, we take the interplanetary dust par- 
ticles (IDPs) as an analog for the dust in debris disks. The 
anhydrous chondritic IDPs collected in the stratosphere possi- 
bly of cometary origin show a highly porous structure (Brown- 
lee 1987). Love et al. (1994) have measured the densities of 

150 unmelted chondritic IDPs with diameters of ^5-15 iim, 
using grain masses determined from an absolute X-ray anal- 
ysis technique with a transmission electron microscope and 
grain volumes determined from scanning electron microscope 
imaging. They found that these particles have an average 
density of ~2.0gcm"^, corresponding to a moderate porosity 
of ~0.4. More recently, Joswiak et al. (2007) identified 12 
porous cometary IDPs (based on their atmospheric entry veloc- 
ities) with an average density of ^ l.Ogcm"^, corresponding to 
P « 0.7. Much higher porosities (>0.9) have been reported for 
some very fluffy IDPs (e.g. MacKinnon et al. 1987, Rietmeijer 
1993), despite that highly porous IDPs are probably too fragile 
to survive atmospheric entry heating. Low densities are also 
derived for different groups of meteoroid^ ascribed to comets: 
the densities are ^ 1.0, 0.75, and 0.27 gcm"^, respectively (Ce- 
plecha et al. 1998), corresponding to a porosity of ~0.7, 0.8, 
and >0.9. 

It is therefore reasonable to assume that the dust in the 
HR4796A disk has a high porosity (at least for those originated 
from cometary bodies). We are not sure about the relative con- 
tributions to the HR4796A disk from asteroid collisions and 
cometary activity. Note that this is a long-standing problem 
even for our own solar system (e.g. see Lisse 2002). 

To conclude, we argue that the presence of tholin-like com- 
plex organic materials in the HR4796A disk is still not con- 
clusive since the observed red spectral shape of the disk can 
be closely reproduced by models of porous dust comprised of 
common cosmic dust species. A more thorough study of the 
scattered light over a range of scattering angles would further 
constrain the optical properties of the dust. 

We thank J.H. Debes and the anonymous referee for their 
very helpful comments. MK and AL are supported in part by 
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AL is supported by the NSFC Outstanding Overseas Young 
Scholarship. 
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Table 1 

Parameters for porous dust models consisting of common dust species (amorphous 
silicate, amorphous carbon, water ice). model 5 is preferred since it closely 

reproduces both the observed reflectance spectrum (see flg.[t]) and the observed ir 
emission (see flg.[3ll. although model 3 fits the observed reflectance spectrum better 
than any other models (see flg.[t]), it does not fit the observed ir emission (see flg-o. 
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Fig. 1 . — Comparison of the model scattered light spectra with the observed visible to near-IR reflectance spectrum of the HR4796A disk of Debes et al. (2008; 
crosses). Diamonds show the model spectra convolved with the HST STIS and NICMOS filters. The model with P = 0.73, a,„i„ = 2 fim and a = 2.9 (red solid line) 
is prefeiTed since it fits both the observed reflectance spectrum and the observed IR emission (see Fig.|3). While the model with P = 0.50, f min = 1 A'lTi and a = 2.8 
(cyan dotted line) fits the reflectance spectrum better than any other models, it does not fit the observed IR emission (see Fig.[3j. In contrast, the model with P = 0.73, 
Qmin = 1 ^^d a = 2.9 (blue dashed line) provides an excellent fit to the observed IR emission, its fit to the observed reflectance spectrum is not as good as the 
other two models shown in this figure. 
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Fig. 2.— Deviations | [fsca(A)/f jj] / [^icaCA)//^;;] - 1 1 of the model results (integrated over flie instrument filters) from the HST STIS and NICMOS 

photometry at each of the seven bands. The deviations are within ~7.5% at all wavebands for the preferred model (with P = 0.73, a,ni„ : 
closely fits both the reflectance spectrum (see Fig.[T) and the IR emission of the HR4796A disk (see Fig.[3). 
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Fig . 3 . — Comparison of the observed IR emission of the HR 4796A dust disk to the model spectra calculated from the porous dust models consisting of amorphous 
silicate, amorphous carbon and ice. 



